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Abstract 

The momentum-weighted sum of the charges of tracks associated to a jet is sensitive to the 
charge of the initiating quark or gluon. This paper presents a measurement of the distri¬ 
bution of momentum-weighted sums, called jet charge, in dijet events using 20.3 fb”^ of 
data recorded with the ATLAS detector at -^ = 8 TeV in pp collisions at the LHC. The 
jet charge distribution is unfolded to remove distortions from detector effects and the result¬ 
ing particle-level distribution is compared with several models. The pr-dependence of the 
jet charge distribution average and standard deviation are compared to predictions obtained 
with several LO and NLO parton distribution functions. The data are also compared to dif¬ 
ferent Monte Carlo simulations of QCD dijet production using various settings of the free 
parameters within these models. The chosen value of the strong coupling constant used to 
calculate gluon radiation is found to have a significant impact on the predicted jet charge. 
There is evidence for a pT-dependence of the jet charge distribution for a given jet flavor. In 
agreement with perturbative QCD predictions, the data show that the average jet charge of 
quark-initiated jets decreases in magnitude as the energy of the jet increases. 
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1 Introduction 


Quarks and gluons produced in high-energy particle collisions hadronize before their electric charge can 
be directly measured. However, information about the electric charge is embedded in the resulting colli¬ 
mated sprays of hadrons known as jets. One jet observable sensitive to the electric charge of quarks and 
gluons is the momentum-weighted charge sum constructed from charged-particle tracks in a jet [1]. Called 
the jet charge, this observable was first used experimentally in deep inelastic scattering studies [2-8] to 
establish a relationship between the quark model and hadrons. Since then, jet charge observables have 
been used in a variety of applications, including tagging the charge of f^-quark jets [9-19] and hadron- 
ically decaying W bosons [20-25] as well as distinguishing hadronically decaying W bosons from jets 
produced in generic quantum chromodynamic (QCD) processes [25] and quark jets from gluon jets [24, 
26-28]. 

The study presented in this paper is a measurement of the jet charge distribution in inclusive dijet events 
from pp collisions at the LHC. Inclusive dijet events provide a useful environment for measuring the jet 
charge as they are an abundant source of gluon-initiated and quark-initiated jets. There are fewer theo¬ 
retical ambiguities from close-by jets and large-angle radiation associated with assigning the jet flavor in 
events with two jets than in events with higher jet multiplicities. Furthermore, the transverse momentum 
(Pt) range accessible in dijet events is broad, (9(10) GeV up to (9(1000) GeV. Since the initial state at 
the LHC has a net positive charge, the probability for positively charged quarks to be produced in pp 
collisions is higher than that for negatively charged quarks. The probability for collisions to involve a 
positively charged valence up-quark in the proton increases with the parton center-of-mass energy V^- 
Thus the average jet charge in inclusive dijet events is expected to increase with V5 if it is correlated 
with the quark charge. The parton distribution functions (PDFs) encode the probabilities to find gluons 
and certain flavors of quarks at given momentum fractions x of the proton. The momentum fractions 
of the two initial partons xi , X 2 , and the proton and parton center-of-mass energies and are re¬ 
lated by = ^Jx\X 2 S. The PDFs are fairly well constrained [29] in the x range relevant for this study, 
0.005—0.5. However, if the jet charge is directly sensitive to the parton flavor, its px dependence can 
provide a consistency check using new information beyond the jet px, which is currently used in PDF 
fits. The PDFs are not the only nonperturbative input needed to model the jet charge distribution and 
its evolution with V5. As a momentum-weighted sum over jet constituents, the jet charge is sensitive 
to the modeling of fragmentation. Previous studies have shown that there are qualitative differences be¬ 
tween the charged-particle track multiplicities of jets in data and as predicted by the leading models of 
hadron production [26]. Thus, a measurement of the jet charge distribution with a range of quark/gluon 
compositions can provide a constraint on models of jet formation. 

While the change in jet parton flavor due to the PDF x-dependence predicts most of the variation in the 
jet charge distribution as a function of V5, there is a second contribution due to the energy-dependence 
of the fragmentation functions. Ratios of the charge distribution moments at different values of V5 
can be calculated perturbatively. Recent calculations [34, 35] in the context of Soft Collinear Effective 
Theory [36-39] show a significant reduction in the magnitude of the average jet charge for a given jet 
flavor as a function of jet energy. Information from PDFs can be used to extract the energy-dependence 
of the average jet charge in the data for direct comparisons to the predictions. 

This paper presents a measurement of the px-dependence of the jet charge distribution’s mean and stan¬ 
dard deviation in dijet events in pp collisions at = 8 TeV with the ATLAS detector. The jet px is a 
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measurable quantity that is strongly related to V^- The average jet charge is extracted for both the lead¬ 
ing and subleading jet and they are distinguished based on their relative orientation in rapidity'. After a 
description of the ATLAS detector (Sec. 2), the data and simulated samples (Sec. 3) and the detector- and 
particle-level objects and selections used in the analysis (Sec. 4), Sec. 5 details the construction of the jet 
charge and some of its properties. In order for the measured jet charge distribution to be compared with 
particle-level models, the data are unfolded to remove distortions from detector effects, as described in 
Sec. 6. Systematic uncertainties in the measured jet charge spectra are discussed in Sec. 7 and the results 
are presented in Sec. 8. 


2 The ATLAS detector 

ATLAS [40] is a general-purpose detector designed to measure the properties of particles produced in 
high-energy pp collisions with nearly a full An coverage in solid angle. The innermost subsystem of 
the detector is a series of tracking devices used to measure charged-particle trajectories bent in a 2 T 
axial field provided by a solenoid whose axis is parallel with the beam direction. This inner detector 
(ID) consists of a silicon pixel detector surrounded by a semiconductor microstrip detector (SCT) and a 
straw-tube tracker. It has full coverage in (p and can detect particles with \r]\ < 2.5. Charged-particle tracks 
are reconstructed from all three ID components, providing measurements of the transverse momentum 
of tracks with a resolution crp^lpT = 0.05% x px/ GeV 0 1%. The track reconstruction algorithm fits 
five track parameters: do, zo, <p, G, and qjp, where do and zo are the transverse and longitudinal impact 
parameters, respectively, q is the track charge and p is the track momentum. Excellent spatial precision 
is required to maintain a well-performing track reconstruction out to and exceeding charged-particle px 
of 1 TeV, where track sagittas are < 0.2 mm. 

Surrounding the ID and solenoid are electromagnetic and hadronic calorimeters to measure showers 
from charged and neutral particles. The high-granularity liquid-argon (LAr) sampling electromagnetic 
calorimeter is located just beyond the solenoid and spans the range Ipl < 3.2. Beyond the electromagnetic 
calorimeter is the two-component hadronic calorimeter that uses scintillator-tile sampling technology in 
the range \q\ <1.7 and LAr sampling technology for 1.5 < \q\ < 3.2. Additional calorimeters are located 
in the forward region. Surrounding the calorimeters is a muon spectrometer, with trigger and precision 
chambers, and incorporating three large toroid magnets composed of eight coils each. 

Due to the large event rate, not every collision can be recorded for processing offline. Events are selected 
using a three-level trigger system [41] that is hardware-based at the first level and software-based for the 
two following levels. An event must satisfy all three trigger levels to be recorded for further processing. 
At each stage of the trigger, energy thresholds are placed on jet-like objects, with the similarity between 
online and offline jets increasing with each level. The first level makes decisions based on low-granularity 
calorimeter towers with thresholds that are typically less than half of the energy required by jets at the 
second level. A simple jet reconstruction algorithm is used at the second level in regions around the 
jets identified by the first level. The third level (known as the Event Eilter) clusters jets with the same 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the center of the detector 
and the z-axis along the beam pipe. The z-axis points from the IP to the center of the LHC ring, and the (/-axis points 
upward. Cylindrical coordinates (r, 0) are used in the transverse plane, (p being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 6 as r/ = -lntan(0/2). The variable AR = -sj^ApY + {ArjY is a 
measure of how close two objects are in the ( 77 , p) plane. The rapidity of a four-vector is defined as 7 / = 0.51n where E 

is the energy and p. the component of the momentum parallel to the beam axis. 
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Trigger threshold [GeV] 

Offline Selection [GeV] 

Euminosity [fb '] 

25 

[50,100] 

7.84x10'^ 

55 

[100, 136] 

4.42x10-4 

80 

[136, 190] 

2.32x10-3 

110 

[190, 200] 

9.81x10-3 

145 

[200, 225] 

3.63x10-2 

180 

[225, 250] 

7.88x10-2 

220 

[250, 300] 

2.61x10-^ 

280 

[300, 400] 

1.16 

360 

>400 

20.3 


Table 1; The single-jet trigger menu used to collect dijet events with the 2012 dataset. The first column is the 
level-three (Event Filter) jet pj threshold and the second column is the offline leading-jet range corresponding 
to the given trigger. The luminosity collected with each trigger is in the last column. The total 2012 dataset was 
20.3 fb“*; the highest-px trigger is not prescaled. 


algorithm as used offline over the entire detector with thresholds that are typically 20-30 GeV higher than 
at level two. The single-jet trigger thresholds increase at each level due not only to differences in the 
jet reconstruction and calibrations, but also to meet the different bandwidth requirements at each trigger 
level. For low-px dijets, the event rate is far too large to save every event that passes the trigger selection 
and so most of the jet triggers are prescaled to artificially lower their recording rate. 


3 Data and simulated samples 

This measurement uses the full dataset o^ pp collisions recorded by the ATLAS detector in 2012, cor¬ 
responding to an integrated luminosity of 20.3 fb“^ at a center-of-mass energy of = 8 TeV. Events 
are only considered if they are collected during stable beam conditions and satisfy all data-quality re¬ 
quirements [42]. To reject noncollision events, there must be a primary vertex reconstructed from at least 
two tracks each with pj > 400 MeV [43] . Due to the high instantaneous luminosity and the large total 
inelastic proton-proton cross section, on average there are about 21 simultaneous (pileup) collisions in 
each bunch crossing. 

A set of single-jet triggers is used to collect dijet events with high efflciency. Table 1 shows the collected 
luminosity for each trigger as well as the offline jet p^ ranges used, chosen such that the trigger is fully 
efficient. The highest-px trigger is not prescaled. 

Monte Carlo (MC) simulated events are generated in pj slices in order to ensure a large number of events 
over a broad range of reconstructed jet px, given constraints on the available computing resources. The 
px slices span the interval 0 to 5 TeV in ranges that approximately double with each increasing slice, 
starting with a range of size 8 GeV and ending with a range of size 2240 GeV. The baseline sample used 
for the measurement is generated with Pythia 8.175 [44] with the AU2 [45] set of tuned parameters (tune) 
and the next-to-leading-order (NLO) PDF set^ CTIO [48, 49]. Another large sample of events is gener¬ 
ated with Herwig-i--i- 2.63 [50, 51] with tune EE3 [52] and leading-order (EG) PDF set CTEQ6E1 [53] 
(particle-level samples with CTIO and EE4 are also used for comparisons). Both Pythia and Hp.RwiG-r-r 


^ A discussion on the use of NLO PDF sets with LO matrix elements is given in Refs. [46, 47]. 
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are LO in perturbative QCD for the (2 —> 2) matrix element and resum the leading logarithms (LL) in the 
parton shower. However, the ordering of emissions in the MC resummation in the shower differs between 
these two generators: Pythia implements px-ordered showers [54] whereas Herwig++ uses angular or¬ 
dering [55]. The phenomenological modeling of the non-pertubative physics also differs between Pythia 
and HFRWTG-t-t. In addition to different underlying-event models (Ref. [56] for Pythia and an eikonal 
model [57] for Herwig-i-i-) the hadronization models differ between Pythia (Lund string model [58]) 
and Herwig-e-i- (cluster model [59]). These two schemes are known [26] to predict different numbers of 
charged particles within jets and different distributions of the charged-particle energies within jets, both 
of which are important for the jet charge. All tunes of the underlying event that are used with Pythia 
and Herwig-e-e in this analysis use LHC data as input. As discussed in Sec. 1, the corrected data are 
compared to models with various PDF sets; for consistency, each set has a dedicated underlying-event 
tune constructed in the same way from a fixed set of data inputs (AU2) described in detail in Ref. [45]. 
The PDF sets include LO sets CTEQ6L1 [53] and MSTW08LO [29] as well as NLO sets CTIO [48, 
49], NNPDF21 NLO [60], and MSTW2008NLO [29]. A sample generated with a NLO matrix element 
from Powheg-Box r2262 [61-64] (henceforth referred to as Powheg) with PDF set CTIO interfaced with 
Pythia 8.175 and the AU2 tune is also used for comparisons. 

All MC samples are processed using the full ATLAS detector simulation [65] based on GEANT4 [66]. 


4 Object reconstruction and event selection 

The reconstructed objects used for the jet charge as well as for the event selection are described in Sec. 4.1. 
The fiducial definition of the measurement, unfolded to particle level, is given in Sec. 4.2. 

4.1 Object reconstruction at detector level 

Jets are clustered using the anti-k^ jet algorithm [67] with radius parameter R = 0.4 implemented in 
EastJet [68] from topological calorimeter-cell clusters [69], calibrated using the local cluster weighting 
(ECW) algorithm [70, 71]. An overall jet energy calibration accounts for residual detector effects as well 
as contributions from pileup [72] in order to make the reconstructed jet energy an unbiased measurement 
of the particle-level jet energy. Jets are required to be central (|? 7 | < 2.1) so that their charged particles are 
within the \t]\ < 2.5 coverage of the ID. 

When more than one primary vertex is reconstructed, the one with the highest ^ of tracks is selected as 
the hard-scatter vertex. Events are further required to have at least two jets with pj > 50 GeV and only the 
leading two jets are considered for the jet charge measurement. To select dijet topologies, the two leading 
jets must have <1.5, where and are the transverse momenta of the jets with the 

highest and second-highest pT^ respectively. The jet with the smaller (larger) absolute pseudorapidity \r]\ 
is classified as the more central (more forward) jet. A measurement of the more forward and more central 
jet charge distributions can exploit the rapidity-dependence of the jet flavor to extract information about 
the jet charge for a particular flavor. This is discussed in more detail in Sec. 4.2. 

Tracks used to calculate the jet charge are required to have px ^ 500 MeV, \t]\ < 2.5, and a per 
degree of freedom (resulting from the track fit) less than 3.0. Additional quality criteria are applied to 
select tracks originating from the collision vertex and reject fake tracks reconstructed from random hits 
in the detector. In particular, tracks must be well-matched to the hard-scatter vertex with |zo sin(0)| <1.5 
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Figure 1; The distribution of the number of tracks associated with a jet in two example jet pj ranges. 


mm and |r/ol < 1 mm, where zo and do are calculated with respect to the primary vertex. Tracks must 
furthermore have at least one hit in the pixel detector and at least six hits in the SCT. The matching 
of tracks with the calorimeter-based jets is performed via the ghost-association technique [73]: the jet 
clustering process is repeated with the addition of ghost versions of measured tracks that have the same 
direction but infinitesimally small pr, so that they do not change the properties of the calorimeter jets. A 
track is associated with a jet if its ghost version is contained in the jet after reclustering. The distribution 
of the number of tracks in two representative jet pr ranges is shown in Fig. 1. The number of tracks 
increases with jet pr and the data fall between the predicted distributions of Pythia and Herwig-i-i-. 


4.2 Object definitions at particle level 

The measurement is carried out within a fiducial volume mafching fhe experimenfal selecfion fo avoid 
exfrapolafion info unmeasured kinematic regions fhaf have addifional model-dependence and relafed un¬ 
certainties. Parficle-level definifions of fhe reconsfrucfed objecfs are chosen fo be as close as possible 
fo fhose described in Sec. 4.1. Particle-level jefs are clusfered from generafed sfable particles wifh a 
mean lifetime t > 30 ps, excluding muons and neufrinos. As wifh fhe defecfor-level jefs, parficle-level 
jefs are clusfered wifh fhe anfi-fc, R = 0.4 algorifhm. In analogy fo fhe ghosf-associafion of fracks fo 
jefs performed af defecfor level, any charged particle clusfered in a parficle-level jef is considered for 
fhe jef charge calculation.^ There musf be af leasf fwo jefs wifh \r]\ <2.1 and pj > 50 GeV. The fwo 
highesf-pT jets must satisfy the same pT-balance requirement between the leading and subleading jet as 
at detector level < 1.5). Due to the high-energy and well-separated nature of the selected 

jets, the hard-scatter quarks and gluons can be cleanly matched to the outgoing jets. While it is possible 
to classify jets as quark- or gluon-initiated beyond leading order in mjet/Sjet [74], the classification is 
algorithm-dependent and unnecessary for the present considerations. In this analysis, the flavor of a jet is 

^ There is no px > 500 MeV threshold applied to charged particles. The impact of applying such a threshold is negligible for 
all px bins except the first two where effects of up to 1% are observed in the mean and standard deviation of the jet charge. 
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Figure 2: For a given jet flavor, (a) shows the fraction / of jets with that flavor in events passing the particle-level 
event selection and (b) shows the pj evolution of the flavor fractions weighted by charge-sign: /up -H /anti-down - 
/anti-up - /down- The forward-central differences between the flavor fractions are largest at low pT, but the highest 
quark-jet purity occurs at high jet p^. The markers for the more forward and central jets are distinguished by their 
blue and red colors, respectively. 


defined as that of the highest energy parton in simulation within a AR < 0.4 cone around the particle-jet 
axis. The jet flavor depends on rapidity and so the two selected jets are classified as either more forward 
or more central; the more forward jet tends to be correlated to the higher-x parton and is less likely to be 
a gluon jet. Figure 2 shows the flavor fraction for the more forward and more central particle-level jets 
passing the event selection. The pj evolution of the sum of the flavor fractions weighted by the sign of 
the parton charge is shown in Fig. 2(b). The forward-central differences between the flavor fractions are 
largest at low pj, but the highest quark-jet purity occurs at high jet pj. 
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(a) (b) (c) 


Figure 3: The (a) the particle-level and (c) detector-level jet charge distribution for various jet flavors in a sample 
of jets with pj > 500 GeV for k - 0.3. (b): the distribution of the jet-by-jet difference between the particle-level 
and detector-level jet charge distributions. The shaded region is used to fit a Gaussian function to extract the bulk 
response resolution, which is cr ~ 0.5 e, where e is the positron charge. 


5 Constructing the jet charge 


There is no unique way to define the jet charge. The most naive construction is to add up the charge of 
all tracks associated with a jet. However, this scheme is very sensitive to lost or extraneous soft radiation. 
Therefore, a weighting scheme is introduced to suppress fluctuations. Using the tracks assigned to a jet 
by ghost association, the jet charge Qj of a jet J is calculated using a transverse-momentum-weighting 
scheme [1]: 


Qj = 


1 

(PtjY 


ieTracks 


( 1 ) 


where Tracks is the set of tracks associated with jet J, qi is the charge (in units of the positron charge) 
of track i with associated transverse momentum pj i, /c is a free regularization parameter, and pjj is the 
transverse momentum of the calorimeter jet. The distributions of Qj for various jet flavors are shown 
in Fig. 3 for /c = 0.3. In the simulation, there is a clear relationship between the jet charge and the 
initiating parton’s charge, as up-quark jets tend to have a larger jet charge than gluon jets. Furthermore, 
gluon jets tend to have a larger jet charge than down-quark jets. However, the jet charge distribution 
is already broad at particle level and the jet charge response (2particie-ievei - Gdetector-ievei) resolution is 
comparable to the differences in the means of the distributions for different flavors, so one can expect 
only small changes in the inclusive jet charge distribution for changes in the jet flavor composition. The 
three narrow distributions on top of the bulk response distribution in Fig. 3(b) are due to cases in which 
only one or two charged particles dominate the jet charge calculation at particle level. The two off-center 
peaks are due to cases in which one of the two high-px-fraction tracks is not reconstructed and the widths 
of the two off-center and central peaks are due to the (single) track and jet px resolutions. The bulk 
response is fit to a Gaussian function with standard deviation cr ~ 0.5 e (units of the positron charge). 

















Figure 4; The distribution of the jet charge built from the leading n tracks {Qin) for (a) k - 0.3 and (b) k - 0.7 in 
a sample of jets with 500 GeV < pi < 600 GeV. The horizontal axis ranges are not the same. In this px range, the 
median number of tracks is about 15. 


The parameter k in Eq. (1) regulates the sensitivity of the jet charge to soft radiation. For /r > 0, the 
jet charge is infrared safe."^ Low values of k enhance the contribution to the jet charge from low-px 
particles while in the k ^ oo limit, only the highest-px track contributes to the sum in Eq. (1). The 
dependence on the highest-px tracks is demonstrated with the plots in Fig. 4 with the variable Qjn, which 
is the jet charge in Eq. (1), but built from the leading n tracks. The variable Qi\ is simply the weighted 
fragmentation function of the leading-track px to the jet px with weight k. The usual Qj is recovered 
in the limit n — > oo. Figure 4 shows the sequence Qjn for k = 0.3 and k = 0.7. For lower values of k, 
many tracks are required for the sequence of distributions to converge to the full jet charge. However, for 
K = 0.7, the distribution converges quickly, indicating that only the highest-px tracks are contributing. 
All reconstructed tracks are henceforth used when computing the jet charge, but the plots in Fig. 4 give 
an indication of the contribution of (relatively) high- and low-px tracks. Dedicated studies [24] agree 
with theoretical predictions [34] that suggest that k ~ 0.5 is the most sensitive to the charge of the parton 
initiating a jet. Therefore, the measurement presented in this paper uses k = 0.5 in addition to /r = 0.3 and 
K = 0.7 in order to maintain a broad sensitivity to both hard and soft radiation inside jets. 

The reconstructed jet charge distributions for k = 0.3 and k = 0.7 are shown in Fig. 5 for events passing 
the selection described in Sec. 4. Section 6 describes how the jet charge moments are corrected for 
detector resolution and acceptance elfects through an unfolding procedure. 


The jet charge is never collinear safe for *: > 0 and not even the sign of the jet charge is Lorentz invariant, although it is clearly 
invariant under longitudinal boosts. 
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Figure 5: The detector-level jet charge distributions for (a,b) 50 GeV < px < 100 GeV and (c,d) 500 GeV < pj < 600 
GeV for the more forward jet and (a,c) k - 0.3 and (b,d) k - 0.7. The peak at zero in the top left plot is due to jets 
without any tracks. 
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6 Unfolding 


The particle-level jet charge distribution’s mean and standard deviation are measured as a function of jet 
Pt- This is accomplished by unfolding a discretized two-dimensional distribution of the jet charge and the 
jet pt and then computing the first two moments of the jet charge distribution in each bin of pj- The jet 
charge distribution is discretized into 15 bins in each of 10 bins of the jet px- The jet charge mean is robust 
against the bias introduced from the discretization procedure, where fewer than five charge bins per px bin 
is required to ensure negligible bias after recovering the mean from the discretized distribution. However, 
the standard deviation of the jet charge distribution is sensitive to the discretization and requires about 15 
charge bins^ in order for the inherent bias due to discretization to be negligible. The jet charge spans the 
range \Qj\ < 1.8 for/c- = 0.3, \Qj\ < 1.2 for a: = 0.5 and \Qj\ < 0.9 for a: = 0.7. Events in the overflow of the 
jet charge distribution are placed in the first or last bins. The px binning is given by: [50,100), [100, 200), 
[200, 300), [300, 400), [400, 500), [500, 600), [600, 800), [800, 1000), [1000, 1200), and [1200, 1500] 
GeV. Figure 6 displays the px-dependence of the jet charge distribution’s mean and standard deviation 
for detector-level data and simulation and for particle-level simulation. The differences between the 
simulated detector- and particle-level distributions give a indication of the corrections required to account 
for detector acceptance and resolution effects in the unfolding procedure. The growing difference between 
the particle- and detector-level average jet charge is due to the loss of charged-particle momentum inside 
jets as a result of track merging. At particle level, the standard deviation of the jet charge distribution 
decreases with increasing px, but at detector level it increases with px due to resolution effects. 

There is no unique way to extract the particle-level distribution of the jet charge and jet px from the 
reconstructed distributions. An iterative Bayesian (IB) technique [75], implemented in the RooUnfold 
framework [76], is used to unfold the two-dimensional jet charge and jet px distribution. In the IB 
unfolding technique, the number of iterations and the prior distribution are the input parameters. In the 
first step, the raw data are corrected using the simulation to account for events that pass the fiducial 
selection af defector level, but not the corresponding selection at particle level; this correction is the. fake 
factor. Then, the IB method iteratively applies Bayes’ theorem using the response matrix to connect the 
prior to posterior at each step, with the nominal Pythia sample used for the initializing prior. The response 
matrix describes the bin migrations between the particle-level and detector-level two-dimensional jet 
charge and jet px distributions. While the response matrix is nearly diagonal, the resolution degrades 
at high Px where more bin-to-bin migrations from particle level to detector level occur. The number of 
iterations in the IB method trades off unfolding bias with statistical fluctuations. An optimal value of four 
iterations is obtained by minimizing the bias when unfolding pseudo-data Herwig-i-i- with Pythia as a 
test of the methodology. The last step of the unfolding applies another correction from simulation to the 
unfolded data to account for the differential rate of events passing the particle-level selection but not the 
detector-level selection; this correction is the inefficiency factor. 


7 Systematic uncertainties 


All stages of the jet charge measurement are sensitive to sources of potential bias. The three stages of the 
measurement are listed below, with an overview of the systematic uncertainties that impact the results at 
each stage: 


^ Fifteen is the number used in the unfolding and differs from the number shown for illustration in e.g. Fig. 5. 
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Figure 6: The detector-level (data and simulation) and particle-level jet charge distribution’s (a) average and (b) 
standard deviation as a function of the jet pj for the more forward jet. The ratios in the bottom panel are constructed 
from the simulation, and show the prediction of detector-level Pythia over the data (top ratio), and detector-level 
Pythia over particle-level Pythia (bottom ratio). Bars on the data markers represent only the statistical uncertainties. 
For both (a) and (b), k — 0.5. 


Correction Factors: Fake and inefficiency factors are derived from simulation to account for the fraction 
of events that pass either the detector-level or particle-level fiducial selection, but not both. These 
factors are generally between 0.9 and 1.0 except in the first pj bin, where threshold effects introduce 
corrections that can be as large as 20%. Experimental uncertainties correlated with the detector- 
level selection acceptance, such as the jet energy scale uncertainty, result in uncertainties in these 
correction factors. An additional source of uncertainty on the correction factors is due to the explicit 
dependence on the particle-level jet charge and jet pj spectra. A comparison of particle-level 
models (Pythia and Herwig-p-p) is used to estimate the impact on the correction factors. 

Response Matrix: For events that pass both the detector-level and particle-level fiducial selections, fhe 
response mafrix describes migrations befween bins when moving befween fhe defecfor level and fhe 
particle level. The response mafrix is faken from simulafion and various experimenfal uncerfainfies 
on fhe jef charge and jef pj specfra resulf in uncerfainfies in fhe mafrix. Uncerfainfies can be divided 
info fwo classes: fhose impacting fhe calorimefer jef pj and fhose impacfing frack reconsfrucfion 
inside jefs. 

Unfolding Procedure: A dafa-driven fechnique is used fo esfimafe fhe pofenfial bias from a given choice 
of prior and number of iferafions in fhe IB mefhod [77]. The particle-level specfrum is reweighfed 
using fhe response mafrix so fhaf fhe simulafed defecfor-level specfrum has significanfly improved 
agreemenf wifh dafa. The modified defecfor-level disfribufion is unfolded wifh fhe nominal re¬ 
sponse mafrix and fhe difference befween fhis and fhe reweighfed parficle-level specfrum is an 
indication of fhe bias due fo fhe unfolding mefhod (in parficular, fhe choice of prior). 
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The following two subsections describe the impact of the detector-related sources of systematic uncer¬ 
tainty in more detail. Uncertainties on the calorimeter jet pT are described in Sec. 7.1 and the uncertainties 
related to tracking are described in Sec. 7.2. Summaries of the systematic uncertainties for the more for¬ 
ward jet and k - 0.5 are found in Table 2 and Table 3 for the average jet charge and the jet charge 
distribution’s standard deviation, respectively.^ The uncertainties for the more central jet are similar. 

Average Jet Charge Jet p^ Range [100 GeV] 


Uncertainty [%] 

[0.5,1] 

[1,2] 

[2,3] 

[3,4] 

[4,5] 

[5,6] 

[6,8] 

[8,10] 

[10,12] 

[12,15] 

Total Jet Energy Scale 

+ 8.4 

- 13.6 

+ 3.8 

- 3.5 

+ 0.9 

- 5.0 

+ 0.8 

- 0.3 

+ 1.1 

- 1.6 

+ 1.1 

- 1.1 

+ 0.7 

- 1.0 

+ 0.7 

- 0.9 

+ 0.4 

- 0.7 

+ 0.9 

- 0.3 

Jet Energy Resolution 

+ 6.8 

- 6.8 

+ 2.3 

- 2.3 

+ 0.7 

- 0.7 

+ 0.7 

- 0.7 

+ 0.3 

- 0.3 

+ 0.3 

- 0.3 

+ 0.1 

- 0.1 

+ 0.1 

- 0.1 

+ 0.1 

- 0.1 

+ 0.3 

- 0.3 

Charged Energy Eoss 

+ 0.0 

- 0.0 

+ 0.0 

- 0.0 

+ 0.0 

- 0.0 

+ 0.0 

- 0.0 

+ 1.7 
- 0.0 

+ 1.5 
- 0.0 

+ 1.5 
- 0.0 

+ 1.5 
- 0.0 

+ 1.6 
- 0.0 

+ 3.6 

- 0.0 

Other Tracking 

+ 3.3 

- 1.6 

+ 0.0 

- 0.4 

+ 0.9 

- 0.2 

+ 0.7 

- 0.1 

+ 0.5 

- 0.4 

+ 1.4 
- 0.6 

+ 0.7 

- 0.9 

+ 1.2 
- 1.2 

+ 1.1 
- 1.3 

+ 0.9 

- 1.7 

Track Multiplicity 

+ 0.0 

- 1.5 

+ 0.1 

- 0.0 

+ 0.0 

- 0.6 

+ 0.0 

- 1.1 

+ 0.0 

- 0.8 

+ 0.0 

- 0.6 

+ 0.0 

- 1.2 

+ 0.0 

- 1.4 

+ 0.0 

- 2.1 

+ 0.0 

- 2.9 

Correction Eactors 

+23 

-23 

+ 0.9 

- 0.9 

+ 0.8 

- 0.8 

+ 1.0 
- 1.0 

+ 0.3 

- 0.3 

+ 0.6 

- 0.6 

+ 0.1 

- 0.1 

+ 0.3 

- 0.3 

+ 0.2 

- 0.2 

+ 0.1 

- 0.1 

Unfolding Procedure 

+28 

-28 

+ 2.4 

- 2.4 

+ 0.3 

- 0.3 

+ 0.2 

- 0.2 

+ 0.2 

- 0.2 

+ 0.3 

- 0.3 

+ 1.1 
- 1.1 

+ 1.0 
- 1.0 

+ 1.6 
- 1.6 

+ 0.6 

- 0.6 

Total Systematic 

+38 

-39 

+ 5.1 

- 4.9 

+ 1.7 
- 5.2 

+ 1.6 
- 1.7 

+ 2.1 

- 1.9 

+ 2.4 

- 1.6 

+ 2.1 

- 2.1 

+ 2.3 

- 2.3 

+ 2.6 

- 3.0 

+ 3.8 

- 3.4 

Data Statistics 

28 

7.4 

1.4 

0.7 

0.3 

0.6 

0.9 

2.0 

4.2 

7.0 

Total Uncertainty 

+47 

-48 

+ 9.0 

- 8.9 

+ 2.2 

- 5.4 

+ 1.8 
- 1.9 

+ 2.1 

- 1.9 

+ 2.5 

- 1.7 

+ 2.3 

- 2.3 

+ 3.0 

- 3.0 

+ 5.0 

- 5.2 

+ 8.0 

- 7.8 

Measured Value [e] 

0.014 

0.024 

0.049 

0.065 

0.076 

0.082 

0.092 

0.100 

0.108 

0.115 


Table 2: A summary of all the systematic uncertainties and their impact on the mean jet charge for k — 0.5 and the 
more forward jet. The correction factors are the fake and inefficiency corrections applied before/after the response 
matrix. The Other Tracking category includes uncertainty on the track reconstruction efficiency, track momentum 
resolution, charge misidentification, and fake track rate. All numbers are given in percent. 


7.1 Calorimeter jet uncertainties 

Jets are calibrated so that the detector-level px is an unbiased measurement of the particle-level jet px 
and various data-driven techniques are used to derive in situ estimates of the difference in this calibration 
between the data and the simulation. Uncertainties in the energy scale and resolution of calibrated jets 
impact the jet charge in the normalization of Eq. (1) (but preserve the jet charge sign) as well as the binning 
for the 2D distribution. Complete details of this source of uncertainty can be found in Ref. [78]. There are 
many components of the jet energy scale uncertainty. The in situ correction is derived from data using the 
momentum balance in events with Z bosons (low px) or photons (moderate px) produced in association 
with jets as well as the balance of multijet (high px) and dijet (high Ipl) systems. Uncertainties on this 
method stem from the modeling of these processes in simulation. There is also a contribution from the 
response to single hadrons [79]. Additional sources of uncertainty are due to the modeling of the in-time 
and out-of-time pileup corrections to the jet energy scale as well as differences in the response due to the 
flavor of the jet. To assess the impact of each component of the jet energy scale uncertainty, the jet energies 

^ The uncertainties on the first pj bin of the average jet charge are much larger than on the other bins because the mean is small 
compared to the resolution. 
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Standard Deviation 

Systematic 
Uncertainty [%] 


Jet pt Range [100 GeV] 


[0.5,1] [1,2] [2,3] [3,4] [4,5] [5,6] [6,8] [8,10] [10,12] [12,15] 


Total Jet Energy Scale 

+ 1.9 
- 1.7 

+ 1.5 
- 1.3 

+ 1.1 

- 1.1 

+ 1.1 

- 1.0 

+ 0.9 

- 0.8 

+ 1.0 
- 0.7 

+ 0.8 

- 0.8 

+ 0.7 

- 0.8 

+ 0.5 

- 0.5 

+ 0.5 

- 0.5 

Jet Energy Resolution 

+ 1.3 
- 1.3 

+ 0.3 

- 0.3 

+ 0.1 

- 0.1 

+ 0.2 

- 0.2 

+ 0.3 

- 0.3 

+ 0.4 

- 0.4 

+ 0.2 

- 0.2 

+ 0.2 

- 0.2 

+ 0.2 

- 0.2 

+ 0.2 

- 0.2 

Charged Energy Eoss 

+ 0.0 

- 0.0 

+ 0.0 

- 0.0 

+ 0.0 

- 0.0 

+ 0.0 

- 0.0 

+ 0.2 

- 0.0 

+ 0.3 

- 0.0 

+ 0.3 

- 0.0 

+ 0.3 

- 0.0 

+ 0.4 

- 0.0 

+ 1.1 
- 0.0 

Other Tracking 

+ 0.0 

- 0.3 

+ 0.1 

- 0.3 

+ 0.2 

- 0.4 

+ 0.3 

- 0.4 

+ 0.4 

- 0.5 

+ 0.5 

- 0.4 

+ 0.5 

- 0.5 

+ 0.5 

- 0.5 

+ 0.5 

- 0.4 

+ 0.4 

- 0.4 

Track Multiplicity 

+ 0.0 

- 0.2 

+ 0.0 

- 0.3 

+ 0.0 

- 0.2 

+ 0.0 

- 0.1 

+ 0.0 

- 0.0 

+ 0.1 

- 0.0 

+ 0.2 

- 0.0 

+ 0.2 

- 0.0 

+ 0.3 

- 0.0 

+ 0.2 

- 0.0 

Correction Eactors 

+ 0.9 

- 0.9 

+ 0.1 

- 0.1 

+ 0.0 

- 0.0 

+ 0.1 

- 0.1 

+ 0.0 

- 0.0 

+ 0.1 

- 0.1 

+ 0.0 

- 0.0 

+ 0.0 

- 0.0 

+ 0.0 

- 0.0 

+ 0.0 

- 0.0 

Unfolding Procedure 

+ 1.9 
- 1.9 

+ 0.4 

- 0.4 

+ 0.0 

- 0.0 

+ 0.1 

- 0.1 

+ 0.2 

- 0.2 

+ 0.0 

- 0.0 

+ 0.1 

- 0.1 

+ 0.3 

- 0.3 

+ 0.4 

- 0.4 

+ 1.7 
- 1.7 

Total Systematic 

+ 3.1 

- 3.0 

+ 1.6 
- 1.5 

+ 1.1 
- 1.2 

+ 1.2 
- 1.1 

+ 1.1 
- 1.0 

+ 1.2 
- 0.9 

+ 1.0 
- 0.9 

+ 1.0 
- 1.0 

+ 1.0 
- 0.8 

+ 2.1 

- 1.8 

Data Statistics 

0.9 

0.3 

0.1 

0.1 

0.0 

0.1 

0.1 

0.3 

0.6 

1.0 

Total Uncertainty 

+ 3.2 

- 3.1 

+ 1.6 
- 1.5 

+ 1.1 

- 1.2 

+ 1.2 
- 1.1 

+ 1.1 

- 1.0 

+ 1.2 
- 0.9 

+ 1.0 
- 1.0 

+ 1.1 

- 1.0 

+ 1.2 
- 1.0 

+ 2.4 

- 2.1 

Measured Value [e] 

0.410 

0.387 

0.375 

0.372 

0.370 

0.369 

0.368 

0.367 

0.362 

0.355 


Table 3: A summary of all the systematic uncertainties and their impact on the jet charge distribution’s standard 
deviation for k - 0.5 and the more forward jet. The correction factors are the fake and inefficiency corrections 
applied before/after the response matrix. The Other Tracking category includes uncertainty on the track reconstmc- 
tion efficiency, track momentum resolution, charge misidentification, and fake track rate. All numbers are given in 
percent. 

in simulation are shifted according to the pt- and Ty-dependent ±lcr variations. For a fixed variation, the 
response matrix, and fake and inefficiency factors are recompufed and fhe unfolding procedure is repealed. 
The resulling uncerfainfy on the jel charge distribution’s mean and standard deviation is about 1% or less 
for jet Pt above 200 GeV. The jet energy resolution uncertainty is derived using data-driven techniques in 
dijet events [80]. To assess the impact of a slightly larger jet energy resolution, jet energies are smeared 
according to px- and p-dependent factors and propagated through the entire unfolding procedure, as for 
the jet energy scale uncertainty. The jet energy resolution uncertainty is subdominant to the jet energy 
scale uncertainty. 

7.2 Tracking uncertainties 

Uncertainties on tracking are broken down into contributions related to (i) the efficiency of reconsfrucling 
Iracks and (ii) measuremenls of Ihose Iracks fhal are successfully reconslructed. The uncertainly on fhe 
inclusive Irack reconslruclion efficiency is dominated by fhe uncertainly in fhe malerial in fhe ID. The 
amounl of material is known to wifhin ~ 5% [81]. Simulated defector geomelries wifh various levels of 
material in fhe ID wifhin fhe measured uncertainties are used fo estimate fhe track reconstruction effi¬ 
ciency uncertainty. These uncertainties are p- and px-dependent, ranging from < 1% for |p| < 2.1 to 
< 4% for 2.1 < IpI < 2.3 and < 7% for 2.3 < |p| < 2.5. The impact of the uncertainty is estimated by 
randomly removing tracks within the px- and p-dependent probabilities, leading to a < 0.5% uncertainty 
on the jet charge distribution’s mean and standard deviation. An additional uncertainty accounts for the 
difference in efficiency befween dala and simulafion due lo fhe modeling of the track per number of 
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degrees of freedom (NDF) requirement. A requirement of;t^^/NDF < 3 is more than 99% efficient across 
jet and track /jt. but the efficiency is generaiiy higher in simuiation than in data. The difference in the 
efficiency between data and simuiation is < i0% of the inefficiency. The impact of this mismodeling is 
evaluated by independently removing tracks with a probability that is 10% of the /NDF < 3 require¬ 
ment inefficiency. As a result of this procedure, the jet charge distribution’s mean and standard deviation 
change by < 0.1% in most pj bins. 

In addition to the loss of tracks due to the material in the ID, tracks can be lost due to the busy envi¬ 
ronment inside the cores of jets. This loss can be studied in simulation by comparing the reconstructed 
charged particle momentum with the charged particle momentum inside the corresponding particle-level 
jet. In order to remove the impact of the tracking resolution and the contribution from fake tracks already 
accounted for separately, reconstructed tracks in the simulation are matched with charged particles. The 
matching is performed by considering the energy deposited in the various layers of the ID by charged 
particles due to material interactions modeled with GEANT4. Weights are assigned to charged particles 
based on the energy deposited in detector elements that were used to reconstruct a given track. A match 
is declared if the weight for one charged particle is sufficiently high. Figure 7(a) shows the ratio in sim¬ 
ulation of the sum of the px of charged particles that were matched to reconstructed tracks to the sum 
of the pt of all the charged particles as a function of the jet px. At low px, the ratio increases with px 
due to losses as a result of hadronic interactions with the material in the ID. Beyond about 200 GeV, the 
fraction monotonically decreases due to the loss of tracks in the core of the jet. A related quantity is 
(E track px / Jet px), where the denominator is the reconstructed calorimeter jet px and the numerator is 
a sum over tracks associated with the jet. Since the particle-level charged-to-neutral fraction of the energy 
is independent of px, a degradation in this ratio can provide information about the loss of tracks inside 
the core of a jet in data. Figure 7(b) shows the distribution of (E track px / Jet px) as a function of jet 
Px. It exhibits trends very similar to those in Fig. 7(a), and in fact the relative loss (fraction with respect 
to the peak) is similar. The MC underestimates the loss by < 1%. The impact of the charged-particle 
momentum loss inside the cores of jets is estimated by randomly removing tracks with a px-dependent 
probability such that the relative loss in the simulation matches that in the data. This uncertainty is negli¬ 
gible for jets with px < 400 GeV, but is non-negligible for higher-px jets, resulting in a < 4% uncertainty 
on the average jet charge in the highest px bin. 

The momentum resolution of isolated tracks has been well measured m JJip pp and Z ^ pp 
events [82]. The scale and resolution of reconstructed muon candidates are shifted and smeared in the 
MC simulation to account for differences between the data and the simulation for Generic tracks 
are not corrected in the same way as muon candidates reconstructed from the ID (and ignoring the muon 
spectrometer). Thus, the correction factors are taken here as the systematic uncertainty on the momentum 
resolution. The momentum resolution is parameterized as a sum in quadrature of a p“^-dependent term, 
a constant term, and a term linear in the track px, with coefficients ro, ri, and r 2 , respectively. The first 
term accounts for fluctuations in the energy loss in the detector material, the second term captures effects 
due to multiple scattering, and the third term accounts for the intrinsic resolution caused by misalignment 
and the finite spatial resolution of ID hits. Unlike muon spectrometer tracks, ID tracks do not traverse a 
significant amount of material and so the energy-loss coefficient vq and its uncertainty are neglected. The 
uncertainties on r \, ri and the momentum scale s are estimated by randomly smearing the px of every 
track with Ipl-dependent factors that are < 2% for track px < 100 GeV and increase to 10-20% at 1 
TeV depending on \r}\. Propagating these variations through the unfolding procedure results in uncertain¬ 
ties that are subdominant to other uncertainties, but non-negligible (~ 2%) in the highest px bins for the 
average jet charge. 
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Figure 7; (a): The average of the px-weighted ratio of charged particles that were matched to reconstructed tracks 
(S Matched charged pj) to all the charged particles (S charged px) as a function of the particle-level jet px- (b): 
The average of the ratio of the sum of px from tracks to the calorimeter jet px (2 track px / Jet px) as a function of 
jet Px in both data and simulation. The momentum ratio of charged particles to all particles is nearly 2/3 due to the 
number of pion species (as indicated by the straight lines for Herwig-h-i- and Pythia predictions at particle level), 
but is not exactly 2/3 due to the presence of photons and kaons in the jet. 


Aside from the track pj, the other track parameter that is relevant for the jet charge is the track charge. 
Especially at high px when the tracks are very straight, the probability for misidentifying the track charge 
can become non-negligible. Simulated particles matched to generator-level information, as described 
above, are used to study the charge flipping probability (charge misidentification rate) for nonfake tracks 
originating from charged particles inside jets (mostly pions). The rate predicted from the simulation is 
<0.1% for track px < 100 GeV, 0.5% for 100 GeV < px < 200 GeV, 1% for 200 GeV < px < 300 
GeV, 2% for 300 GeV < px < 400 GeV and 4% for px > 400 GeV. Dedicated studies of charge flipping 
in searches for same-sign leptons [83] suggest that the mismodeling of the charge-flipping rate is (much) 
less than 50%. Therefore, the impact of charge flipping on the jet charge measurement is conservatively 
estimated by randomly flipping the charge of tracks at 50% of the charge misidentification rate. The 
impact on the measured jet charge mean and standard deviation is negligible. 

Random combinations of hits in the detector can be combined to form a reconstructed track. Tracks 
resulting in particular from multiparticle trajectories that have kinks can result in very large reconstructed 
track Px- The quality criteria are effective at mitigating the presence of fake tracks, which constitute less 
than < 0.1% of all reconstructed tracks. To determine the impact of fake tracks on the jet charge, fake 
tracks are randomly removed with a probability that is 50% of the rate in simulation. This results in a 
negligibly small uncertainty on the mean jet charge and a < 0.5% uncertainty on the standard deviation 
of the jet charge distribution. 

The tracking uncertainties described so far take into account the resolution and efficiency of reconstruction 
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of charged-particle momenta. One last source of systematic uncertainty is the number of charged particles. 
The unfolding procedure uncertainty takes into account the uncertainty on the prior due to the charged- 
particle multiplicity, but the jet charge resolution also changes with the charged-particle multiplicity. To 
assess the impact on the response matrix of the mismodeled charged-particle multiplicity, the distribution 
of ntrack is reweighted in the simulation per pj bin and the relative dilference when unfolding the nominal 
Pythia distribution with the reweighted Pythia distribution is taken as a systematic uncertainty.^ This 
uncertainty is subdominant for the standard deviation across pj and for the mean at low to moderate jet 
Pf. For the mean jet charge, the largest uncertainty is with the smallest ic and for large pi, where it is 
3-4% percent in the highest px bin for k = 0.3 and k = 0.5. 


8 Results 

The data satisfying the event selection criteria described in Sec. 4 are unfolded according to the procedure 
in Sec. 6 and the average and standard deviation of the jet charge distribution are computed as a function 
of the jet px- These results, along with the systematic uncertainties detailed in Sec. 7, are discussed in 
Sec. 8.1. The PDF uncertainty and jet formation uncertainties in the theory predictions are compared to 
the unfolded data in Secs. 8.2 and 8.3, respectively. Using PDF information as input, the average charge 
per jet flavor is extracted in Sec. 8.5 and its px-dependence is studied in Sec. 8.6. 

8.1 Unfolded Jet Charge Spectrum 

The unfolded jet charge mean is shown as a function of the jet px in the top plots of Fig. 8 for /c = 0.3, 
0.5 and 0.7. The average charge increases with jet px due to the increase in up-flavor jets from PDF 
effects. The average charge increases from 0.01c at px ~ 100 GeV to 0.15c at px ~ 1.5 TeV. Systematic 
uncertainties are generally a few percent, except at low jet px where the fractional uncertainty is large 
because the average jet charge in the denominator is small, and at high px where the tracking uncertainties 
are not negligible. The statistical uncertainty is estimated by repeating the measurement on an ensemble 
of bootstrapped datasets: each event is used in each pseudodataset n times, where n is a random number 
distributed according to a Poission distribution with mean one. The first bin suffers from large statistical 
uncertainties (up to 170%), but for the higher px bins the systematic uncertainty is dominant, except at 
the highest px bin where statistical and systematic uncertainty are of similar size (about 7%). The jet 
charge distributions of the more forward and more central jet differ in shape, in particular at low px, due 
to the different shape of the up/down flavor fractions in those bins as shown in Fig. 2(b). 

Analogous results for the standard deviation of the jet charge distribution are shown in the bottom plots 
of Fig. 8. Even though the standard deviation of the reconstructed jet charge distribution increases with 
jet Px (Fig. 6), the particle-level value decreases and approaches an asymptote for px > 300 GeV. 


8.2 Sensitivity of PDF Modeling 

Variations in the PDF set impact the relative flavor fractions and thus in turn change the jet charge distri¬ 
bution. Such changes do not vary much with k, since the PDF impacts the jet charge distribution mostly 


’ Since the prior is also changed, this uncertainty at least partially includes the unfolding procedure uncertainty. 
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Figure 8: The measured average (standard deviation) of the jet charge distribution on top (bottom) in units of the 
positron charge as a function of the jet px for k - 03,0.5, and 0.7 for the more forward jet (left) and the more 
central jet (right). The crossed lines in the bars on the data indicate the systematic uncertainty and the full extent of 
the bars is the sum in quadrature of the statistical and systematic uncertainties. 
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Figure 9: The average jet charge (k - 0.5) in units of the positron charge for (a) the more forward jet and (b) the 
more central jet compared with theory predictions due to various PDF sets. The crossed lines in the bars on the 
data indicate the statistical uncertainty and the full extent of the bars is the sum in quadrature of the statistical and 
systematic uncertainties. 


through the flavor fractions. Figures 9 and 10 compare the unfolded distributions of the jet charge distri¬ 
bution’s average and standard deviation with several PDF sets, with tuned predictions for Pythia for each 
PDF, and with the same AU2 family of tunes. The sampling of PDF sets results in a significant spread for 
the average jet charge, but has almost no effect on the standard deviation. CTEQ6L1 describes the data 
best, although the data/MC ratio has a stronger pj dependence. In particular, the data/MC differences 
with CTEQ6L1 are up to 10% (15%) at moderate pT for the more forward (central) jet. Eor high pj, 
differences between data and simulation are less significant. NEO PDEs such as CTIO are consistently 
below the data by about 10%-15%. 


8.3 Sensitivity of QCD Models and Tunes 

The measurements presented in Sec. 8.1 show that there are qualitative differences between the data and 
the MC simulations, and comparisons in Sec. 8.2 suggest that variations in the PDE set cannot fully ex¬ 
plain the differences. Differences in Sec. 8.1 between Pythia and Herwig-f-f suggest that some aspect of 
the modeling of fragmentation could lead to the observed differences between the simulation and the data. 
One possible source is the hadronization modeling, which differs between Pythia (Eund-string fragmen¬ 
tation) and Herwig-f-f (cluster fragmentation). The modeling of final-state radiation (ESR) is expected to 
have an impact on the jet charge distribution because variations in the radiation lead to different energy 
flow around the initial parton and hence different fragmentation of the jet. The plots in Eig. 11 and Eig. 12 
show the measured average jet charge and the jet charge distribution’s standard deviation, respectively, 
for K = 0.3, 0.5, and 0.7, compared to various models for a fixed PDE set (CTEQ6E1). In addition to 
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Figure 10: The standard deviation of the jet charge (k - 0.5) distribution in units of the positron charge for (a) 
the more forward jet and (b) the more central jet compared with theory predictions due to various PDF sets. The 
crossed lines in the bars on the data indicate the statistical uncertainty and the full extent of the bars is the sum in 
quadrature of the statistical and systematic uncertainties. 


Pythia 8 and Herwig++ model predictions, Figs. 11 and 12 contain the predictions from Pythia 6 using 
the Perugia 2012 tune [84] and the radHi and radLo Perugia 2012 tune variations. These Perugia tune 
variations test the sensitivity to higher/lower amounts of initial- and final-state radiation (via the scaling 
of as), although only variations of the FSR are important for the jet charge distribution. For the mean jet 
charge, Pythia 6 with the P2012 radLo tune is very similar to Pythia 8 with the AU2 tune. The spread 
in the average jet charge due to the difference befween fhe radHi and radLo funes increases wifh k, since 
suppression of soff radiation makes fhe jef charge disfribufion more sensifive fo fhe modeling of fhe en¬ 
ergy fraclion of fhe leading emissions. For fhe jef charge disfribufion’s sfandard deviation, fhe sensifivify 
fo fhe as scaling is large af bofh high and low k. However, fhe sensifivify is inverfed: radHi gives a larger 
sfandard deviation for k = 0.3, buf a lower sfandard deviafion for k = 0.7. Ofher Perugia 2012 funes 
have been sfudied, fesfing fhe sensifivify fo color-reconnecfion and multiple parfon inferacfions, buf fhe 
differences in fhe jef charge disfribufion’s mean and sfandard deviafion are small. The Perugia 2012 funes 
may nof fully capfure fhe spread in nonperfurbafive effecls, which is also suggesfed by fhe increasing 
difference befween Pythia 8 and Herwigj-j- for decreasing k. 
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Figure 11: The average of the jet charge distribution in units of the positron charge for (a) k - 0.3, (b) 0.5, and (c) 
0.7 comparing various QCD MC models and tunes for the more forward jet. The crossed lines in the bars on the 
data indicate the statistical uncertainty and the full extent of the bars is the sum in quadrature of the statistical and 
systematic uncertainties. 
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Figure 12: The standard deviation of the jet charge distribution in units of the positron charge for (a) k - 0.3, (b) 
0.5, and (c) 0.7 comparing various QCD MC models and tunes for the more forward Jet. The crossed lines in the 
bars on the data indicate the statistical uncertainty and the full extent of the bars is the sum in quadrature of the 
statistical and systematic uncertainties. 
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8.4 Model comparison overview 


Figures 13 and 14 show comparisons of the unfolded jet charge distribution’s mean and standard deviation 
for different QCD simulations using LO and NLO PDF sets. The predictions using the CTIO NLO PDF 
set as shown in Fig. 13 are generally about 10% below the data. Consistent with the expectation that 
the PDF and (nearly collinear) fragmentation are responsible for the jet charge distribution’s mean and 
standard deviation, there does not seem to be an effect from the Powheg NLO matrix element. For the 
jet charge distribution’s standard deviation and k - 0.3, the data falls between PYTHIA (larger standard 
deviation) and Herwig++ (smaller standard deviation), but this trend is less evident for larger k values, 
suggesting a difference due to soft tracks. As seen in Sec. 8.2, comparisons with CTEQ6L1 show it be to 
a better model for the pT-dependence of the mean jet charge than CTIO. The analogous plots to Fig. 8 
but using CTEQ6L1 instead of CTIO are shown in Eig. 14. Generally, there is agreement between the 
simulation and the data with only a < 5% difference in the lower pT bins. 


23 


0 . 8 - 


hiJl 


*♦* 



i*’ jtt 


500 




1000 1500 

Jet p [GeV] 



0 . 8 - 




tt? »» #f i»‘ 

500 


_lU 

1000 1500 

Jet [GeV] 


fO 

■> 

o 

Q 


0.8 


0.6 


0.4 


0.2 


ATLAS 

- is = 8 TeV, 20.3 fb 

- More Forward Jet 

- PDF: CTIO 

Data 

— Pythia8.175AU2 _ 

- - Herwig++ 2.63 EE4 

Powheg + Pythia 8.175 AU2 


* K = 0.3 • K = 0.5 ^k = 0.7 




- 







cc 

■> 

0 

Q 



1.1 

1 

0.9 

1.1 

1 

0.9, 



** * A A. 



: 


Tj *1 > S N ^ 

s 

s 

; 


Jet p_ [GeV] 


hi*’.*'*•’*•’*** *•’ ^ 

500 1000 


1500 

Jet [GeV] 


Figure 13: The measured average of the jet charge distribution (top), and the standard deviation (bottom), in units 
of the positron charge as a function of the jet pj for k - 0.3,0.5, and 0.7 for the more forward jet (left) and the more 
central jet (right) using CTIO as the PDF set. The markers in the lower panel are artificially displaced horizontally 
to make distinguishing the three k values easier. The Powheg+Pythia curves are nearly on top of the Pythia curves. 
The crossed lines in the bars on the data indicate the systematic uncertainty and the full extent of the bars is the sum 
in quadrature of the statistical and systematic uncertainties. 
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Figure 14: The measured average of the jet charge distribution (top), and the standard deviation (bottom), in units 
of the positron charge as a function of the jet pj for k — 03,0.5, and 0.7 for the more forward jet (left) and the 
more central jet (right) using CTEQ6L1 as the PDF set. The markers in the lower panel are artificially displaced 
horizontally to make distinguishing the three k values easier. The crossed lines in the bars on the data indicate 
the systematic uncertainty and the full extent of the bars is the sum in quadrature of the statistical and systematic 
uncertainties. 
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8.5 Extraction of the average up-quark and down-quark jet charges 

In addition to understanding the trends in the jet charge distribution from PDFs, one can use PDFs to 
extract information about jets of a particular flavor. These exclusive interpretations rely on flavor-fraction 
information in PDFs and matrix element calculations to extract the jet charge distribution for particular 
jet (anti-)flavors in each px bin. The required nonperturbative information is summarized in Fig. 2(a). 
Jets with flavors other than up/down/anti-up/anti-down/gluon are not included in Fig. 2(a) and give a 
negligible contribution (< 2%) in the highest px bins. 

One way of extracting the up- and down-flavor average jet charges is to exploit the difference in flavor 
fractions shown in Fig. 2(a) between the more forward and the more central jets. Due to the px-balance 
requirement between the leading and subleading jet in the event selection, to a good approximation, the 
px spectrum is the same for the more forward and the more central jet. Assuming that the average jet 
charge of the sum of flavors that are not up/down/anti-up/anti-down is zero, in each bin i of px: 


/ ^forward \ _ / ^forward _ ^forward'\ i /^forward _ ^forward x^down 

~ \Jup,i anti-up,i/^-^down,! anti-down,i'^^i 

/ ^central \ _ / ^central _ ^central \ ^'^P i /^central _ ^central x^down 
'^7 ~ Vup,! anti-up,i/^•^down,i •^anti-down,i^^i ’ 

where Qj is the jet charge from Eq. (1), f^. is the fraction of flavor y in px bin i for the jet x € 
{more forward, more central} and Q^. is the average jet charge for such jets. The values f^. are taken 
from simulation (Pythia with CTIO PDF and AU2 tune), which then allows an extraction of by solv¬ 
ing the system of equations in Eq. (2). This extraction is performed separately in each px bin. Eigure 15 
shows the extracted up- and down-flavor jet charges in bins of jet px. At very high jet px, the absolute 
quark flavor fractions are large (Eig. 2), but the difference between the more forward and more central jets 
is small and the statistical uncertainty is large. At low jet px, the difference between the more forward 
and more central jets is large (Eig. 2), but the absolute quark flavor fraction is small and the statistical 
uncertainty is once again large because the mean jet charge is close to zero. In the limit that the flavor 
fractions are identical for the more forward and more central jet, the equations become degenerate and it 
is not possible to simultaneously extract the average up- and down-flavor jet charges. The uncertainties 
on the flavor fractions and on the measured average jet charges are propagated through the solutions of 
Eq. (2). Generally, the uncertainty is larger for the down-flavor jets because the fraction of these jets is 
smaller than the fraction of up-flavor jets. 


8.6 Dependence of the up-quark and down-quark jet charge on px 

Using the methods of Sec. 8.5, one can examine the residual px-dependence of the average jet charge 
after accounting for PDE effects. The inclusive jet charge has been shown to increase with px due to a 
mixing of jet flavors and the following subsection investigates the px-dependence of a fixed jet flavor. 
Section 8.6.1 describes the theory predictions and the extracted px-dependence from the data is discussed 
in Sec. 8.6.2. 
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Figure 15: The extracted value of up- and down-quark jet charges in units of the positron charge in bins of jet pi 
for K - 03,0.5, and 0.7. The error bars include statistical, experimental systematic, and CTIO PDF uncertainties 
added in quadrature. The thick part of the error bar indicates the PDF contribution to the total uncertainty and the 
horizontal line on each error bar indicates the contribution from the statistical uncertainty. The first two pi bins are 
excluded due to their very large uncertainties. 
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8.6.1 Theory prediction 


Recent theoretical studies have shown that the energy dependence of jet charge moments is calculable 
perturbatively [34, 35]. At leading power (treating /njet/f'jet as the expansion parameter), the charge de¬ 
fined with jet and track energies is equivalent to the one given in Eq. (1), since the jet opening angle is 
small. For jets defined by a radius R af an energy E, fhe average jef charge (Qj) is given by 


(Qj) = [1 + 0(as)] 2] QhD’^iK, E X R), (3) 

h 

where Qu is fhe charge of hadron h and fhe funcfions D^{k, fi) are fhe Mellin momenfs of fhe fragmenfafion 
functions 





dx^D^^^{x,n). 


(4) 


The fragmenfafion funcfions /r) describe fhe probabilify for a hadron h fo carry a momenfum fraction 
X of a quark q af fhe energy scale [85, 86]. The (9(as) correction is small and is dominafed by fhe 
uncerfainfy on fhe fragmenfafion funcfions [34]. Ratios of {Qj) af differenf energies resulf in a cancellation 
of fhe leading corrections and associafed uncerfainfies. Soff corrections fo fhe average jef charge are small 
nof only because of fhe addifional suppression af /c > 0 compared fo collinear radiafion, buf also because 
fhe leading soff emissions are made of gluons, which carry no elecfric charge. Alfhough gluons can splif 
info quark and anfi-quark pairs, fhe same number of quarks and anfiquarks go info fhe jef so fhaf on 
average fhe jef charge is unchanged.^ 

The leading energy dependence of fhe average jef charge is due fo fhe derivafive of fhe fragmenfafion 
funcfions wifh respecf fo E, which is defermined by fhe renormalizafion group equafions for D^(A',/r): 


+ 0{al), (5) 

where Pqq{K) is fhe momenf of fhe leading-order spliffing funcfion: 

— 1 + 

PqqiK) = Cf I dziz" - 1)--. (6) 

Jo - z 

Recalling fhaf E/pt is consfanf af leading power for all parficles in fhe jef, fhe predicfion for fhe scale 
violation of a quark jef is given by 


Pt d 

{Qk) dpi 


(Q^) = —PqqiK) = 

7T 


-0.024 ± 0.004 K = 0.3 
■ -0.038 ± 0.006 K = 0.5 
-0.049 ± 0.008 K = 0.7 


(V) 


This is not strictly true, as the soft radiation pattern depends on the partonic colors, which are correlated with electric charge, 
but these correlations are negligible compared to the other uncertainties. 
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where the last form consists of numerical approximations setting ji = E x R equal to 50 GeV and 500 
GeV with their average giving the central value, and using cTsCSO GeV) = 0.130 and as(500 GeV) = 
0.094. Experimentally, one measures combinations of quark and gluon jets with the fractions of different 
partonic flavors varying with energy. The next subsection discusses how the scaling violation parameter, 
defined in Eq. (7), can be extracted from the data. 


8.6.2 Extraction from the data 

Since from Eq. (7), one can approximate a linear dependence on c^: 


{Qj){pi) = <2(1 + G ln(pT/PT)) + 0{cl), (8) 

where Q - {Qj){p^) for some fixed (but arbitrary) transverse momentum, px- Therefore, for a fixed pj 
bin i, the measured charge is given as a superposition of the average jet charge for various jet flavors: 


(Qi) « + Gln(pT.,/PT)), (9) 

/ 

where Pfj is the fraction of flavor / in bin i, 2/ is the average jet charge of flavor / and pj is a fixed 
transverse momentum. Eitting the model in Eq. (9) directly to the data to extract Qf is not practical 
because there are three parameters and only 10 px bins, some of which have very little sensitivity due 
to low fractions /? or large uncertainties on (Qj). One way around this is to extract Qf in one fixed bin 
of transverse momentum (denoted px) as described in Sec. 8.5. Then Eq. (9) is highly constrained, with 
only one parameter for which each other bin of px gives an estimate. The systematic uncertainties are 
propagated through the fit treated as fully correlated between bins and the statistical uncertainty is treated 
coherently by bootstrapping^. A weighted average is performed across all px bins and for both the more 
forward and the more central jet. The procedure is summarized below: 

1. In the bin 600 GeV < px < 800 GeV, extract the values 2up and 2down- These values can be seen 
in the fifth px bin of Eig. 15. 

2. With 2up and 2down fixed, extract the scale violation parameter estimate c*.,; in each px bin i by 
solving 

(2i)measured “ + c^j ln(px,;/px)) (10) 

./■ 

where px = 700 GeV is the bin center from the previous step. 

3. Repeat the above procedure for all systematic variations and for all bootstrap pseudo-datasets to 
arrive at estimates of the uncertainty (T{cK,i) for each px bin i. 


^ Pseudo-datasets are generated by adding each event in the nominal dataset j times where j is a Poisson random variable with 
mean 1. Since events are coherently added, this respects the correlations in the statistical uncertainty for the more forward 
and central jet charges. 
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Figure 16: The extracted values of the scale violation parameter q from the data compared to theoretical calcula¬ 
tions [34, 35]. The error bars include statistical, experimental systematic, and PDF uncertainties added in quadra¬ 
ture. The thick part of the error bar indicates the PDF contribution to the total uncertainty and the horizontal line on 
each error bar indicates the contribution from the statistical uncertainty (each shown without adding in quadrature 
any other source of uncertainty). 


4. The central value for the extracted scale violation parameter is - (X; CKjl(T{Ciij))l li,(l/o‘(c;f,)). 

5. The uncertainty cr{cK) is determined by repeating step (3) with the nominal values replaced by 
their systematic varied versions or the bootstrap pseudo-data values for the statistical uncertainty 
estimate. 

The results are presented in Fig. 16. The data support the prediction that < 0 and dc^ldK < 0. Linear 
correlations between k values can be determined using the bootstrapped datasets: about 0.9 between co .3 
and C0.5 as well as between cq.j and cq.v, while the correlation is about 0.7 between co .3 and cqj. Thus, 
the three points are quite correlated, but there is additional information from considering more than one k 
value. 
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9 Summary 


This paper presents a measurement of the particle-level px-dependence of the jet charge distribution’s 
mean and standard deviation in dijet events from 20.3 fb“^ of = 8 TeV pp collision data recorded 
by the ATLAS detector at the LHC. The measured jet charge distribution is unfolded to correct for the 
detector acceptance and resolution for direct comparison to particle-level models. Comparisons are made 
at particle level between the measured jet charge distribution and various PDF sets and models of jet 
formation. Simulations with Pythia 8 using the CTEQ6L1 PDF set describe the average jet charge of 
the more forward jet within about 5% and the more central jet within about 10%. The jet charge dis¬ 
tribution’s standard deviation is described within 2%. HuRwio-r-r shows a similarly good agreement for 
Pt > 500 GeV and k = 0.7. However, the Herwig-i-i- predictions decrease systematically for both the 
average and the standard deviation for decreasing k. Predictions with the CTIO NLO PDF are systemat¬ 
ically below the data across jet px for the average jet charge and systematically above for the jet charge 
distribution’s standard deviation. Taking the PDFs as inputs, the average up- and down-flavor jet charges 
are extracted as a function of px and are compared with predictions for scale violation. The data show 
that the average up- and down-quark jet charges decrease slightly with px and this decrease increases 
with K, as predicted. The particle-level spectra are publicly available [87] for further interpretation and 
can serve as a benchmark for future measurements of the evolution of nonperturbative jet observables to 
validate QCD MC predictions and tune their free model parameters. 
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E. Conventi'"^®’^, M. Cooke'*, B.D. Cooper**, A.M. Cooper-Sarkar'*", T. Cornelissen'**, M. Corradi*"®, 
E. Corriveau**’’^, A. Corso-Radu'"*, A. Cortes-Gonzalez'*, G. Cortiana'"', G. Costa"'®, M.J. Costa'"*, 
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MJ. Da Cunha Sargedas De Sousa^^^‘‘’^^^'’, C. Da Via^^, W. Dabrowski^^^, A. Dafinca^^°, T. Dai^^, 

O. Dale^^, F. Dallaire^^, C. Dallapiccola^^, M. Dam^^, J.R. Dandoy^^, N.P. Dang^*, A.C. Daniells^^, 

M. Danninger^^*, M. Dano Hoffmann^^^, V. Dao"^^, G. Darbo^°^, S. Darmora^, J. Dassoulas^, 
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I. Deigaard'*’^, J. Del Peso*^, T. Del Prete^^^‘‘’^^^’’, D. Delgove^^’, E Deliot^^^, C.M. Delitzsch^^, 

M. Deliyergiyev^^, A. Dell’Acqua^*^, E. Dell’Asta^^, M. Dell’Orso^^^®’*^^'’, M. Della Pietra^*’^®’-^, 

D. della Volpe^^, M. Delmastro^, PA. Delsart^^, C. Deluca^°^, D.A. DeMarco^^^, S. Demers^’^, 

M. Demichev^^, A. Demilly^*’, S.P Denisov^^*’, D. Derendarz^^, J.E. Derkaoui^^^‘*, E. Derue^°, 

P Dervan^"^, K. Desch^', C. Deterre^^, PO. Deviveiros^*^, A. Dewhurst^^\ S. Dhaliwal^^, 

A. Di Ciaccio^^^^*’^^^*’, E. Di Ciaccio^, A. Di DomenicoC. Di Donato^°^‘*’'°^’’, A. Di Girolamo^®, 

B. Di Girolamo^*^, A. Di Mattia^^^, B. Di Micco^^^^’^^^'’, R. Di Nardo^^, A. Di Simone^^, R. Di Sipio^^^, 

D. Di Valentino^^, C. Diaconu^^, M. Diamond^^^, EA. Dias"^^, M.A. Diaz^^^, E.B. Diehl^^, J. Dietrich^^, 

S. Diglio*^, A. Dimitrievska^^, J. Dingfelder^^, P Dita^^'’, S. Dita^^'’, E. Dittus^®, E Djama*^, 
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T. Dohmae^^^, J. Dolejsi^^^, Z. Dolezal^^^, B.A. Dolgoshein^^’*, M. Donadelli^^‘^, S. Donati^^^^’^^^*’, 
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E. Eilthaut^^^ M. Eincke-Keeleri^'^, K.D. Einelli‘5°, M.C.N. Eiolhais^2^'‘’^2®^ E. Eiorini^^^ A. Eiran^^, 
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E. Eullana Torregrosa*^, B.G. Eulsom^^^, T. Eusayasu^°^, J. Euster^^^, C. Gabaldon^^, O. Gabizon'^^, 
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M.P Geisler^^®, C. Gemme^°^, M.H. Genest^^, S. Gentile^^^^’^^^'’, M. George^^, S. George^^, 
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J. Goncalves Pinto Firmino Da Costa^^^, F. Gonella^^, S. Gonzalez de la Hoz^^^, G. Gonzalez Parra^^, 
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V. S. Jabbar^, B. Jackson^^^, M. Jackson^^, P. Jackson^ M.R. Jaekel^*^, V. Jain^, K. Jakobs^^, 
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A. Karamaoun^, N. Karastathis^*’’^®^, M.J. Kareem^^, E. Karentzos^*^, M. Karnevskiy^^, S.N. Karpov^^, 
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P Klimek'^^®’^^^*’, A. Klimentov^^, R. Klingenberg^^, J.A. Klinger^^^, T. Klioutehnikova^*’, 

E.-E. Kluge^^®, P Kluit^®’, S. Kluth^‘’\ J. Knapik^^, E. Kneringer^^, E.B.EG. Knoops^^, A. Knue^^, 

A. Kobayashi^^^, D. Kobayashi^^^, T. Kobayashi^^^, M. Kobel"^^, M. Koeian^"^^, P Kodys^^^, T. Koffas^^, 
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R. Kowalewski^^^, T.Z. Kowalski^^®, W. Kozaneeki^^^, A.S. Kozhin^^°, V.A. Kramarenko^^, 
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